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We investigate a diffusive motion of a system of interacting Brownian particles in quasi-one- 
dimensional micropores. In particular, we consider a semi-infinite ID geometry with a partially 
^Sl ■ absorbing boundary and the hard-core inter-particle interaction. Due to the absorbing boundary 

the number of particles in the pore gradually decreases. We present the exact analytical solution 
f^ ' of the problem. Our procedure merely requires the knowledge of the corresponding single-particle 

^Sl , problem. First, we calculate the simultaneous probability density of having still a definite number 

» , ■ N — k oi surviving particles at definite coordinates. Focusing on an arbitrary tagged particle, we 

C^ ' derive the exact probability density of its coordinate. Secondly, we present a complete probabilistic 

description of the emerging escape process. The survival probabilities for the individual particles 
are calculated, the first and the second moments of the exit times are discussed. Generally speaking, 
although the original inter-particle interaction possesses a point-like character, it induces entropic 
repulsive forces which, e.g., push the leftmost (rightmost) particle towards (opposite) the absorbing 
boundary thereby accelerating (decelerating) its escape. More importantly, as compared to the 
reference problem for the non-interacting particles, the interaction changes the dynamical exponents 
which characterize the long-time asymptotic dynamics. Interesting new insights emerge after we 
interpret our model in terms of a) diffusion of a single particle in a A'^-dimensional space, and b) 
order statistics defined on a system of A'^ independent, identically distributed random variables. 
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I. INTRODUCTION 

Diffusion dynamics of hard-core interacting particles in one dimension (ID) is known as the single-file diffusion 
(SFD). The SFD occurs for instance in narrow biological pores[l[, inside the zeolites0, y|, during the sliding of 
proteins along the DNA[J, Q. Other examples comprise migration of adsorbed molecules on surfacesQ, diffusion in 
nanotubes[7|,|8|, carrier migration in polymers and superionic conductors [9|, [lOj . Advances in technology now allow to 
study artificial single-file systems such as solid-state nanopores|ll]. confined colloid particles |13 - [l4| . charged spheres 
in circular channel[ia|, or millimetric steel balls in plane condenser [Ifl]. In all these systems, the mutual passage of 
the particles is forbidden which results in deviation of tagged particle's dynamics from the classical diffusion. 

The concept of SFD was first introduced in the mathematical literature by Harris [ITJ. He proved that the mean- 
square displacement of the tagged particle increases with time as i^ " in contrast to its linear increase for the free 
particle. For a comprehensible review of the present state of the art in the theory of SFD cf. Introduction in Ref. 
\l&l ) . Despite the long history of SFD phenomena, exactly solvable models of single-file systems are rare and they have 
occurred only recently. The exact solution of the diffusion equation for N hard-core interacting particles has been 
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2C). The SFD of iV particles 



the present authors have 



first derived in Ref. [19| for the diffusion on an infinite line (it has been rederived in Ref. 
within the finite interval has been solved through the Bethe ansatz in Ref. [2l|. In Ref. 
generalized the previously known exact solutions to the case of an arbitrary time- and space-dependent force acting on 
particles. All these exact solutions can describe real physical situations, however, they all concern the systems with a 
conserved number of particles. Often the opposite is true and the single-file systems should be treated as open systems 
with variable number of particles. Many such systems are listed in the first paragraph. The single-file conditions 
in combination with a variable number of particles are two prerequisites for a proper understanding of kinetics of 
diffusion-limited catalytic reaction in crowded environments 23j . Also the search for corrections to Smoluchowski 
theory which would incorporate the effects of excluded volume |24l [25| stimulates the interest in one-dimensional 
models l26| with hard-core interaction and particle absorption. 

To our knowledge, the only work that presents analytical study of dynamics of a tagged particle in an open system 
is Ref. |l9| where the finite interval with two perfectly absorbing boundaries was briefly discussed. However, the 
approximative tagged particle probability density function (PDF) obtained in the thermodynamic limit was only 
considered. In the present paper we develop a new general formalism that allows an exact analytical treatment of the 
SFD problem with the absorbing boundary (such boundary can represent, e.g., the open end of the biological pore, a 
fixed reactant, a binding site inside the pore, etc.). The approach is a logical outgrowth of our earlier method[22l| and 
its characteristic feature is that the exact solution of the SFD problem can be given whenever the exact solution of the 
corresponding single-particle problem is known. Formally, our analysis is based on a hierarchical system of coupled 
generalized diffusion equations describing the evolution periods when the pore contains a definite number of particles. 
We were able to solve the hierarchy. The solution yields the complete information concerning the many-particle 
dynamics. The physical consequences of the particle interaction manifest themselves most transparently if we follow 
the motion of an individual tagged particle. We were able to derive the exact expression for the tagged particle PDF 
and for its survival probability. The hard-core inter-particle interaction induces entropic forces which radically alter 
the character of the escape process comparing to the corresponding problem of noninteracting particles. 

The paper is organized as follows. Sec. HIl contains the definition of the model. Sec. IIIII comprises four subsections. 
First, for the sake of a latter reference, Subsec. IIII Al presents the complete solution of the underlying single-particle 
case. In Subsec. IIIIBI we formulate the hierarchy of coupled generalized diffusion equations for the partial PDFs. 
The appropriate choice of the initial conditions is discussed in Subsec. IIII CI In Subsec. IIIIDI the SFD dynamics is 
translated into the diffusion of a single compound particle in a special iV-dimensional domain. Sec. IIVI contains the 
exact solution of our model (some purely technical steps are particularized in appendixes). The Eqs. (j24p and (1251) in 
Subsec. IIV Al give the exact solution of the hierarchy of coupled diffusion equations from which the tagged particle's 
PDF (Eq. ([55)1 ) is derived in Subsec. IIVBI We then turn, in Sec. |Vl to a detailed analysis of the escape process. 
Subsec. FV Al discusses the survival probabilities and their asymptotic behavior. The exit times are analyzed in Subsec. 

II. DEFINITION OF THE MODEL 

Consider N identical Brownian particles, each with the diffusion constant Z?, diffusing in the one-dimensional 
semi-infinite pore. The particles are subject to the following two interactions: 

(i) The particle-particle hard-core interaction. As the consequence, particles cannot pass each other. Their initial 
ordering xn < xn-i < ... < Xi is preserved during the subsequent diffusion process. 

(ii) The particle-boundary interaction. We assume a partially absorbing boundary at a; = 0. The particle that hits 
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FIG. 1. An illustration of the one-dimensional semi-infinite pore containing N = G diffusing particles and the partially absorbing 
boundary at the origin (the exit from the pore). The diffusing particles interact as hard spheres and hence their initial ordering 
xe < X5 < ... < a::2 < a;i is retained during the time evolution. 

the boundary is either removed from the pore, or it is immediately reflected and it continues to diffuse inside the 
pore. The boundary will be implemented through the radiation boundary condition [27]. The "particle-capture 
capacity" of the boundary is characterized by the parameter k. The value k = corresponds to a reflecting 
boundary, the limit k — )■ cxa yields a totally absorbing boundary. 

The particles cannot enter the pore from the outside and they can leave it only through the boundary at the origin. 
The number of particles inside the pore either decreases or remains constant over time. 

III. EVOLUTION EQUATIONS 

A. Single-particle case 

In order to introduce basic concepts needed in the following, let us first consider the simplest situation iV = 1. 
That is, at the initial time t — 0, there is only one particle in the pore and it is located at the position y, y > 0. The 
initial PDF p{x; 0\y;0) for the particle position is the Dirac delta function 

p{x;0\y;0) = S{x-y). (1) 

Inside the pore (i.e., for a; > 0), the time-evolution of the initial condition ([T|) is governed by the diffusion equation 

^^p{x;t\y;0) = D-^p{x;t\y;0). (2) 



At the boundary, we require the radiation boundary condition (BC) 

= . (3) 
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The exact solution of Eq. ^ with the initial condition ([1]) and BC ([3]) reads 
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where erfc(») stands for the complementary error function [28J. It can be derived, e.g., by the Laplace transform 
methodflil. 

Let /(x), a; > 0, be an arbitrary PDF describing the initial particle position. Having calculated the Green function 
(j4|) , the time evolution of the initial condition is given as 

p[x;t)^ / dyp(a;;i|2/;0)/(y) . (4) 

Jo 

The probability of finding the particle at the time t within the interval (0, x) is given by the cumulative distribution 
function 

F{x,t) = I Ax'p(x';t) . (5) 

Jo 

The probability that the particle has not exited the pore by the time t, i.e., its survival probability, reads 

S{i) = lim F{x,t) . (6) 

a;— ^+00 

It is an important feature of our setting that all subsequent A^-particle results can be formulated through the single- 
particle functions p(x;t), F{x,t), S(t), and only through them. 



A'^-particle case 



Suppose that, at the initial time t — 0, there are N particles in the pore and that their positions < xn < xn-i < 
... < xi, cf. Fig. 1, are distributed according to the joint PDF /(^'(a?^;0), if at = (xi, ... ,xjv)- Having in mind that the 
total number of particles inside the pore is not conserved, we introduce partial (iV — fc)-particle PDFs f^^~^'{xiq-k]t), 
k = 1,... ,N — 1, xjv-fc = {xi,...,XN-k)- We give them the following meaning: f^^~''^XN^k',t)dxN~k is the prob- 
ability that, at the time t, there are N — k particles in the pore (this is stressed by the superscript) and that these 
particles are located in [xM-k^XN-k + dxN-k)- Altogether, we have N partial PDFs. Further, let /^^•'(t) be the 
probability that, at the time t, there are no particles left in the pore. The knowledge of the partial PDFs and of the 
probability f^^' (i) represents the complete probabilistic description of the resulting single-file diffusion problem. 

The partial PDF f^^^{xN',t) varies in time due to the following two reasons: a) the motion of N interacting 
particles, b) the possible exit of the particle No. N out of the pore. The diffusion equation 






together with the "non-crossing" boundary conditions 
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(7) 
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i = 1,...N — 1, account for the diffusion of N hard-core interacting particles inside the pore. The non-crossing 
boundary conditions ([5]) can be derived by introducing the A^-particle probability current and by assuming that 
the normal component of the current vanishes at the hyperplanes {xn : < xn < ... < Xn = Xn-i < •■■ < xi < -t-oo}, 
n = 2,...,N. Their derivation and the further discussion can be found in Ref. |23. Eqs. (3)-(5). BCs ([5]) ensure that 
the initial ordering of the particles is preserved. Hence, it is only the leftmost particle (the particle No. N) that can 
leave the pore. The exit process is incorporated by the radiation BC (cf. Eq. ([3])) 
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The diffusion equation ([7|) must be supplemented by an initial condition f^^'{xN',0). This will be done in Subsec. 
llHCl 

Let us now focus on the partial PDF f'-^~^\xi\[-i',t). It varies due to the following three reasons: a) the motion 
of A^— 1 interacting particles, b) a possible exit of the particle No. A^ — 1 out of the pore, c) the possible exit of the 
particle No. A^ while the remaining particles are situated in the infinitesimal domain {xn^i, xn-i +dxAr_i). We thus 
arrive at an inhomogeneous generalized diffusion equation with a source term on its right-hand side: 






dt 



(10) 



The last term on the right-hand side accounts for the above point c). The hard-core interaction of A^ — 1 particles is 
again incorporated through non-crossing BCs 
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j — 1, ... N ~ 2. The radiation BC describing the possible escape of the particle No. N —1 reads 
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(11) 
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Presently, the initial condition supplementing Eq. (|10|) is simply f'^^~'^\xN-i]Q) = 0. This follows from the proba- 
bilistic interpretation of PDF f'^^~^^XN-i]t) and from the assumption that, at t = 0, there are A^ particles inside 
the pore. 

Following a similar reasoning, the generalized diffusion equation that governs the evolution of the partial PDF 
/(^-'=)(fjv-fc; i). A: = 2, 3, . . . , A^ - 1, reads 



^ f^^^-'HxN-klt) ^dY: ^f^^-'H^N-k-.t) + n /(^-'=+i)(f^^,+i;i) 



dt 
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Again, the equation is supplemented by the BCs 

/ d d 

\dxj dxj+i 

j — 1, ... N — k — 1, by the radiation BC 
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and by the initial condition f^^~'^\xM-k\ 0) — 0. Of course, in the case k — N ~ 1^ there is no non-crossing BC. 

In the last step of the recursive scheme we focus on the time evolution of the probability /*•"-' (t). It increases due 
to the escape of the last remaining particle and it obeys the simple evolution equation 



dt 
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(16) 



with the initial condition /'•^-'(O) — 0. 

Summing up the above steps, we have obtained a hierarchy of N + 1 coupled evolution equations ([7]), (I13p for 
A: = l,..., A^— 1, and ([T6l). The diffusion equation (O is closed. Together with Eqs. (|8]), ^ and with the initial 
condition f^^'{xN',0) it constitutes a well-defined initial-boundary value problem with a unique solution, the partial 
PDF f^^'{xN;t). Suppose that we know it. Then we insert its value at xnI^ ^q ~ (■^i' ■•■' ^n-i,0) into the source 
term on the RHS of the inhomogeneous generalized diffusion equation (|10[) . Then the inhomogeneous equation (|10l) . 
supplemented by Eqs. ((TT|) . (rT2]) . f'^^~^^XN-i',0) = constitutes again a closed initial-boundary value problem. 
Having obtained its solution, the partial PDF /'^^~^^(xAr_i; i), we again evaluate it at afjv-ilj, ^g and use as the 

source term in the equation for partial PDF f^^~'^\xN^2] t). The hierarchy is terminated by the Eq. (|16p . 



C. Initial distribution of particles 



Among all possible forms of initial distributions there exists one which leads to a considerable simplification of 
the subsequent analytical calculations. Even more importantly, the distribution is physically quite natural for real 
single-file systems. It describes the state which emerges, e.g., after previous constitution of a steady state ^], or after 
previous autonomous relaxation towards a thermodynamic equilibrium. Namely, in this paper we consider the initial 
joint PDF for the A^-particle system in the form 
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for xm ^ TIn , 



(17) 



where TZn denotes a phase space defined as a space of all possible configurations of A^ hard-core interacting particles 
on the semi-infinite interval [0, -t-oo), i.e.. 



TZn 



= {xn ■ < xpf < xn-1 < ■■■ < xi < 



-oo} 



(18) 



The symbol xn G Ti-N means that the components of the vector xn — (xi, ... , xm) respect all the inequalities in (jlSp . 
If xn ^ Ti-N, then at least one of them is violated. 

In general, the function f{x) in Eq. (fT7|) can be any PDF defined and normalized on the half-line x e [0, 4-oo). The 
normalization of the function f{x) implies the proper normalization of the initial condition (J17l) . This can be checked 
by applying Eq. (jA2[) from Appendix |X] where we have summarized the basic rules concerning the integration over 
the phase space TZjy. 

In the following, we derive also the marginal PDF for the position of the rt-th particle. At the time t = 0, it is 
obtained by integrating the function (TTTI) over the coordinates of all other A^ — 1 particles. If ?/;„ (x„; 0) denotes the 
marginal PDF in question, we have 



V'r)(x„;0)= / [dfA,\a;„]/W(fA,;0) 
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where Nl — N ~ n {Nn = n — 1) stands for the number of particles focated left (right) of the n-th particle. The 
symbol [dx^Xxn] is explained in Eq. (|A3p from Appendix 1X1 and the symbol TZ^Xxn, designating the integration 
domain, is explained in the paragraph just above Eq. ()A3| . 

The last expression has been written in the form which suggests two important interpretations. First, suppose 
that we have generated N points on the half-line, each of them being independently drawn from the distribution 
with PDF f{x). Then ipn (x; 0) is the PDF for the n-th rightmost of them (i.e., for the n-th order statistics 30]). 
Second, consider the system of N identical, indistinguishable, and noninteracting particles. Let the position of each 
of them be described by PDF f{x). Then T^n \x; 0)da; equals the probability that there is a particle in (x, x -f Ax) 
and, simultaneously, there are precisely N^ {Nr) particles at arbitrary coordinates to the left (to the right) of it. 

In the following, provided we need a further specification of PDF /(a;), we always take 

/(a;) = ^exp(-|) , L>Q. (20) 

This expression together with Eq. P7)) yields a complete specification of the initial state. The state can be prepared 
in the following way. Up to the time t = 0, a) the pore is closed (reflecting boundary at the left end of the pore 
at a; = 0), b) each particle is pushed towards the reflecting boundary by the external space-homogeneous force (e.g. 
the particles are charged and we assume an external electric field), c) the iV-particle system has reached the state of 
thermal equilibrium. At the time t = instantaneously (adiabatically) a) the pore becomes opened, i.e., we replace 
the reflecting boundary by the partially absorbing one, b) the external field is switched off. 

D. Mapping on single- particle diffusion in N dimensions 

The dynamics described in the previous subsection possess a noteworthy geometric interpretation. It resides in 
mapping of the N interacting diffusing particles onto a single "compound" particle in N dimensions. At t = there 
are N particles inside the pore. Until the exit of the A^-th particle from the pore, the situation corresponds to the 
diffusion of the compound particle in the N dimensional wedge domain 

{xn : < xat < xn-1 < ■■■ < X2 < a;i < +00} 

bounded by A^ — 1 reflecting hyperplanes 

{xn '■ < Xn < ... < Xn = Xn-1 < ■■■ < Xi < +(X)} , 

n = 2, ..., A, and by the absorbing hyperplane 

{xn ■ Xn —0,0 < a^TV-i < ■•• < a;2 < xi < +00} . 

The diffusion within this wedge is described by the partial PDF f'^^'>{xN]t) and it is governed by the A^-dimcnsional 
diffusion equation ([7]). In this picture, the non-crossing boundary conditions ([8]) describe the reflection of the compound 
particle on A^ — 1 reflecting hyperplanes. They guarantee the zero probability currents in the directions perpendicular 
to these hyperplanes. The exit of the left-most particle from the pore corresponds to "sticking" of the compound 
particle to the absorbing hyperplane and it is incorporated by the requirement ([9|). 

The exit of the A^-th particle translates to the further diffusion of the compound particle on the above absorbing 
hyperplane. Equivalcntly, it diffuses within the A— 1-dimcnsional wedge domain 

{xn-1 ■ < Xn-1 < ■■■ < X2 < a;i < H-oo} 

bounded by A — 2 reflecting hyperplanes 

{xn-1 ■ < Xn-1 < ... < Xn — 2:„_i < ... < xi < +00} , 

n = 2,...,A— 1, and by the absorbing hyperplane 

{xn-1 ■ Xn-1 = 0,0 < xn-2 < ■■■ < 2:2 < xi < -foo} . 

The diffusion within this domain is described by the partial PDF f'-^^^^{xN-i',t), etc. 

Eventually, after the exit of the A — 2 particles, there are two last particles left inside the pore. This case is 
illustrated in Fig. [2] where the compound particle diffuses in 2D wedge domain {x2 : < 0:2 < xi < -l-oo} bounded 
by just one reflecting boundary {x2 : < 0:2 = xi < -l-oo} and by the absorbing half-line {x2 : 0:2 = 0, < xi < -t-00}. 
Its motion is described by the partial PDF f^'^\x2;t). After this event, the diffusion of the compound particle is 
described by partial PDF f^^^ {xi; t). The exit of the left particle corresponds to the sticking of the compound particle 
to the half-line {x2 : 0:2 = 0, < 2:1 < +00}. The exit of the last particle simply means the sticking of the compound 
particle to the origin. The probability of finding the compound particle at the origin at t is given by the function 
f'^'^^it). For a further discussion of the 2D case cf. Fig. [51 
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FIG. 2. Mapping of the two-particle SFD onto the diffusion of one compound particle in the 2D wedge domain. The upper 
panel shows the space-time representation of the diffusion within the pore. The red line a; = is the absorbing boundary. 
Whenever the particles collide, their trajectories are mutually reflected. One such event occurs roughly at the time t « 1.6. 
The second particle leaves the pore approximately at the time t ~ 3.2. At the moment, the coordinate of the first particle 
is roughly xi ~ 4.1. The lower panel presents the diffusion of the compound particle within the wedge domain bounded by 
the reflecting boundary (the blue line) and by the absorbing boundary (the red line). The point A= [1,0.5] shows the initial 
position of the compound particle. Whenever its trajectory hits the blue line it is reflected. The reflection corresponds to the 
collisions of the two particles. At the point [a:;i,a;2] ~ [4.1,0], the trajectory of the compound particle sticks to the absorbing 
boundary. This represents the exit of the second particle from the pore at the time t « 3.2. Afterwards, the compound particle 
continues to diffuse along the line X2 = 0. Its final position at the time t = 3.5 is shown by the point B« [5.5, 0]. 



IV. EXACT SOLUTION OF EVOLUTION EQUATIONS 



Partial PDFs 



In our recent work we have shown that the Green function (GF) for the diffusion equation of N hard-core interacting 
particles diffusing in an arbitrary external potential can be constructed from the GF p{x] t \ y; t') of the corresponding 
single-particle problem (see Eq. (6) in Ref. |22). In the present context, the same construction yields GF that solves 
Eq. ([7]) with the non-crossing boundary condition ([5]), with the radiation boundary condition ([H]), and with the initial 
condition p--^'{xN', t' \ ijN', t') — 5{xn — Vn)- The explicit form of this GF is given by the permanent ISlj 



P^'^XxN-MvN^t') 
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if xjv, Vn € Ti-N^ and it vanishes otherwise. The summation runs over all A'^! permutations of initial positions. Having 
the GF, the initial condition ([T7|) evolves as 



/W(xAr;i) = /dy^pW(f^;i|yjv;0)/W(yjv;0) . 



(21) 



TIn 



In the present paper we always assume the initial condition p7|) and it is obvious that this special form simplifies the 
analytical expressions. As a matter of fact, using PT|) . the above function reduces to the simple product (cf. Appendix 



f'^^HxN;t) ^ Nlp{xN;t)p{xN~i;t)...p{xi;t) , 



(22) 



where p{x; t) is the known solution of the corresponding single-particle problem (cf. Eq. @). Notice that the product 
yields the solution just within the domain TZn- It turns out that the products of this type will play an important role 
in the following. Hence we define 
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Nl JJ p(x„; t) for xn E TZn 
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for Xn ^ TZn 



We are now prepared to formulate the exact solution of our diffusion problem. 

Suppose that the initial PDF f^^'{xN;0) possess the form (J17p . Then the hierarchy of coupled evolution equations 
as formulated in Subsec. Illl Bl vossesses the solution 



f(^-'\xN.,;t) = {^^]p(^~'^HxN^,;t) 



nf dt'p{0;t') 





for k^O,l,...,N- 1, and 



/(«)(t) 



t -,N 

kI dt'p{0;t') . 
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Turning to the proof, the solution of the first element of the hierarchy (the case fc = 0) has already been derived 
above. Actually, provided there are still N particles in the pore, the solution ((M)) coincides with Eq. (1^^ . 

For an arbitrary fixed k = 1,. . . ,N — 1, we again first consider the (iV— fc)-particle GF that solves Ec^. p^ with 
the non-crossing BC (ITil) . with the radiation BC (|15l) . and with the initial condition p'^^~''\xN~k',t' \yN-k',t') — 
5{xN~k - VN^k)- It reads 



(N-k)\ 



p^^ ''\xN-k\t\yN-k;t')^ ^ p{xi]t\y^^^C^y,t')p{x2\t\y^^(2)\'t')-P{xN-k;t\y^^(N-ky,i'), 



(26) 



if XN^k,yN-k G Ti-N-k, and it vanishes otherwise. 

In order to derive the partial PDF f^^~''\xN-k,t) we first use the above GF (|26l) to evolve the initial condition 
/(^~'''\xjv-fc; 0). This yields the first term on the right-hand side of the following equation. Secondly, we take into 
account the source term in Eq. p^ which gives the second term on the right-hand side of the following equation. 
Altogether, we obtain 



/(^-'=\fjv.fc;t)= J dijN-kP^''-''\xN-k;t\m-k;0)f^''-''\m-k;0) + 

Am-k I At' p^^-'^) {xN-k ; 1 1 yN-k-, t') K /(^-'=+i) (yjv-fe+i ; 
/ Jo 



(27) 



t') 



VN^k+l=0 



The first term on the right-hand side vanishes because f'^^~''\xj^_k;0) — (initially, there are N particles in the 
pore). As for the second term, we start with fc = l and we invoke the same reasoning which enables to rewrite the Eq. 
(I^TI) in the simple product form (P^ . More explicitly, inserting the product form ([^^ with xjv = into Eq. I^E^ and 



repeating the steps as described in Appendix [B] we carry out the required integration over the space variables. After 
these steps we get 



f^''-'\xN-i;t) ^ Np^^~'\xN^i;t)K dt'p{0;t') 
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/(^-2)(fAr_2;t) = NiN - l)p^''-^\xN-2;t)K [ dt'p{0;t') K f dt"p{0;t") 

Jo Jo 

which equals (cf. Eqs. ([MIIMI) ) 



t -|2 

kI dt'p{0;t') 



(28) 



This is just Eq. ([M)) for fc = 1. The procedure can be simply repeated for fc = 2. After introducing Eq. ([^5]) into Eq. 
we have 



(29) 



This proves the statement (p4)) for fc= 2 The proof is completed by repeating the procedure for fc = 3, . . . , iV. 



(30) 



B. Tagged particle's PDF 

In the preceding section, we have derived the probabilities f^^~^\xN~k',t)dxN-k of the compound event that, at 
the time t, the pore contains N—k particles and that they are found in {xN-k,XN-k + dxN-k)- We now deduce the 
properties of the single, say n-th particle (we shall call it the n-th tagged particle) . We start from the partial PDFs 
/(^"'^^(zTV-fc; t), we integrate out the coordinates of all other particles and we sum over all fc. 

The probability that, at the time i, the pore contains N—k particles and the n-th particle is found in the infinitesimal 
vicinity of the coordinate a;„ is i/'n ^ {xn ', t) da;„ , where (for the definition of the symbols used cf. the text after Eq. 
p9)) and the last paragraph in the Appendix \Q 



V'i^-'^-)(x„;t) =y [dx;N-k\Xn]f^''-''\xN-k;t) 



(31) 



To find the n-th particle in the pore, the pore must contain at least n particles. Hence the marginal PDF of the n-th 
particle's position regardless the positions of all other particles and regardless the number of the particles in the pore 
is given by 



N- 



Mx;t)^ ^^i^-'=)(x;i). 



(32) 



fe=0 



We now insert the partial PDFs as given by ([M)) into Eq. (I5T]) . Skipping out a considerable amount of purely algebraic 
steps, the PDF for the n-th tagged particle emerges in the form 



ipnix;t) = 



N\ 



(n- l)!(iV-n)! 



p{x;t) 



+ 00 



dx'p{x';t) 



n-l 



[1-S{t)+F{x,t)] 



N-r. 



(33) 



Notice that, as promised above, the result is formulated solely by means of the single-particle functions p{x;t), S{t), 
and F{x,t), as defined in Eqs. (jll), ([5]), and dU, respectively. In closing the calculation we would like to emphasize 
the relation 



N 



Y^Mx;t) ^ Np{x;t) 



(34) 



which can proved by the direct summation of the expressions (j33p . The result is also plausible if we apply the reasoning 
elaborated in the following paragraph. 

In the text following Eq. ([T^ we have emphasized the analogies of our SFD model with a) the system of N 
noninteracting indistinguishable particles, and with b) the distribution of the n-th order statistics based on the 
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FIG. 3. The survival probabilities S„{t) for the tagged particles (Eq. (|39|) . solid blue lines), and the single-particle survival 
probability S{t) (Eq. ([SSp . the dashed red line) as the functions of time. The parameters used: N = Q, L = 5m, D = lm^s~^. 



1 ms . The number n ■ 



i on the right near the blue lines denotes the order of the corresponding particle, cf. Fig. 1. 



sample of N independent, identically distributed random variables. However, in Subsec. IIII C[ the discussion refers 
to the initial distribution of the particles. Eq. p3l) testifies that the analogies are valid for an arbitrary time if we 
only supplement them with the necessary modifications stemming from the escape process. Indeed, ipn{x;t)dx gives 
the simultaneous probability that one of the N non- interacting indistinguishable particles is found in {x,x + dx), 
n — 1 particles are located in (x, +oo), and each of the remaining N — n particles is either already removed from the 
pore (with the probability 1 — S{t)) or it is still diffusing in the domain in (0,x) (with the probability F{x,t)). The 
indistinguishability is behind the factorials. They account for all possible sequences of particles to the right and to 
the left of the tagged particle. On the other hand, consider N independent, identically distributed random variables 
each one possessing the probabihty density p{x-, t) + (1 — S{t)) S{x). Then Eq. ([55]) is just the probability density for 
the 71-th biggest among them. In this second analogy, the time enters as a parameter. 



ESCAPE PROCESS 



A. Survival probabilities 



If we use the initial condition ()20|) and if we carry out the steps as given in Subsec. lilIA| the single-particle survival 
probability ([6]) reads 



Sit) = 



1 



kL- D 



(Dt\ JVDt 
kL exp I — I ertc — - — 



— D exp K 



D 



crfc 




(35) 



At t = 0, the survival probability equals to one. It approaches zero as i — >■ oo, i.e., the particle will certainly exit the 
pore. The survival probability psp decays with time as the power law|32l| 



S{t)^ At-^ 



as t ^ oo 



(36) 



This asymptotic time dependence is determined solely by the geometry of the problem (ID semi-infinite interval with 
the absorbing boundary). The constant factor A includes the parameters L,D,k. It reads 



1 



L 



D 



(37) 



We shall demonstrate that the parameters L, D, k enter the time-asymptotic behavior of the interacting particles 
only through this combination. 

Having prepared the single-particle survival probability S(t) we now turn to the system of interacting particles. 
The survival probability Sn (t) of the n-th tagged particle is defined as 



Sn{t) ^ dx^pn{x;t) , n=l,2,...,N , 
Jo 



(38) 
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where ipn{x;t) is the n-th particle's PDF as given by Eq. ((33|) . The main steps of the emerging calculation are 
summarized in Appendix [C] The resulting expression includes only the single-particle survival probability: 

This is the main formula of the present section. 

To get a better insight into the escape process in question, let us examine the difference between the survival 
probabilities of the adjacent particles. We have 

Sn{t) - 5„+i(t) = (^^) {sityr (1 - s{t)f-^ . (40) 

The difference is simply the probability that the particle No. n + 1 has already exited the pore and, simultaneously, the 
particle No. n is still diffusing inside the pore. At the same time, the right-hand side by itself suggests an interpretation 
in terms of a system of noninteractin q a nd indistinguishable particles (the survival properties of N noninteracting 
particles are discussed e.g. in Refs. l33l - l35l ). Eq. (PH)) tells us that the following two events have the same probability: 
(i) exactly TV — n of the initial number N of interacting particles have already exited the pore by the time t, and (ii) 
exactly A^ — n of the initial number N of noninteracting and indistinguishable particles have already exited the pore 
by the time t. Differently speaking, the mass transport out of the pore is not affected by the hard-core interaction. In 
particular, in the probabilistic sense, it will take the same time until all N particles escape from the pore regardless 
they interact or not. 

Fig. 13] illustrates the tagged particle survival probability Sn{t) for different n = 1,...,6, together with the single- 
particle survival probability S{t). Several observations are worth emphasizing. First of all, it follows from Eq. (j40p 
that the function 5„(t) are ordered as 

Si{t) > S2{t) > ... > SN{t) , (41) 

where the equalities hold at t = (all A^ particles are initially inside the pore, thus all S'„(0) = 1) and in the limit 
t — > oo. Each hard-core particle acts on the adjacent ones as a moving (fluctuating) reflecting boundary. For example, 
the exit of the rightmost particle is hindered by the remaining A^ — 1 particles which greatly increases its survival 
probability Si{t). Effectively, the collisions between the particles induce entropic repulsive forces. The leftmost 
particle is pushed by the remaining ones towards the pore's boundary which significantly reduces its lifetime inside 
the pore; its survival probability S'Ar(i) rapidly decreases to zero. As a direct consequence of Eq. (p4)) we obtain 



1 ^ 
-Y^Sr.{t) = S{t) (42) 



which relates the arithmetic mean of tagged particles' survival probabilities to the survival probability of the single- 
diffusing particle. Hence Si{t) always exceeds S{t) but their difference Si{t) — S{t) vanishes as i — >■ co. 

The above consideration manifest themselves most pronouncedly if we consider the time-asymptotic behavior of 
the survival probabilities. Each S'„(t) tends to zero, i.e., each particle will certainly exit the pore. The long-time 
asymptotics of the individual Sn (t) can be obtained directly from Eq. (pQ)) . We have 

Sn{t) - f^jA'H-^ , as t^c5o , (43) 



where A has the same meaning as in Eq. (j36|) . In words, each survival probabilities exhibits again a power-law 
decay. However, the decay exponent depends on the particle's order n and it does not depend on the total number of 
particles. The parameter A^ enters just through the constant prefactor of the power law. This can be understood on 
physical grounds. Notice that the survival probability of the rightmost particle Si{t) decays as i~^/^, i.e., with the 
same dynamical exponent as S{t) (cf. Eq. ((36|) '). Indeed, in the long-time limit it is highly probable that all other 
particles have already escaped and thus the initially rightmost particle behaves as a free one. This is also why the 
total number of particles A^ does not enter the exponents of S'n(t), for an arbitrary n. After a "long enough" time, the 
N—n particles which were initially to the left of the n-th one have already escaped. Thus, irrespective of A^ , the n-th 
particle behaves as if it were the leftmost one in the system of n particles. The total number of particles A^ specifies 
only the expression "long enough" in the above reasoning. 

Finally, let us return to the mapping which has been developed in Subsec. flll Dl where we have discussed the analogy 
between our SFD problem and the diffusion of a single compound particle in the A^-dimensional space. Initially, the 
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compound particle departs from a general interior point in the A''-dimensional wedge domain. The exit of the leftmost 
particle has been translated to the adsorption of the compound particle onto the (A'^— l)-dimensional absorbtion 
hyperplane. Then the compound particle continues to diffuse within this hyperplane until its motion is further 
restricted onto the (A^ — 2)-dimensional absorbtion hyperplane, and so on. The function S'„(i) equals the probability 
that, at the time t, the motion of the compound particle has not been restricted onto the (n— l)-dimensional hyperplane, 
yet. That is, it still moves in some fc-dimensional wedge, where k > n. Similarly, the difference S'„(i) — Sn+i{t) as 
given in Eq. (j40[) yields the probability that, at the time t, the compound particle diffuses in the interior of the 
n-dimensional wedge domain. On the whole, as a byproduct, we have solved a nontrivial problem concerning the 
geometrically restricted A^-dimensional diffusion. 

B. Exit times 

The survival probability Sn{t) (Eq. ((39|) ) provides complete information about the escape process of the n-th 
particle. The PDF 0„(f) of the n-th particle's exit time is given bv|27| 

dSnjt) 
(Pnit) = ^^ . (44) 

Hence all the moments of the exit time can be obtained from S'n(t). It follows from the long-time properties of Sn(t) 
(cf. Eq. (gSl)) that 

Mt) -'^(^)A"t-^-\ as t^oo. (45) 

Consequently, the mean exit time of the n-th particle, 

(■-l-oo 

(t„)= / dtt(j)n{t) , (46) 

is infinite for n < 2. The integral (|46| becomes finite for n > 2 and its convergence does not depend on the total 
number of particles initially located in the pore. TV only modifies the magnitude of the mean exit time (if the latter 
exists). Similarly, the second moment of the exit time (i^) only exists for n > 3. 

The above analysis brings us to an interesting comparison between the escape process confined by the static and by 
the fluctuating reflecting boundaries, respectively. Considering the semi-inflnite interval with the absorbing boundary 
at the origin, the survival probability S{t) of the single-diffusing particle decays to zero as the power law ((36|) and 
the mean exit time is infinite. If we add the static refiecting boundary on the right of the single-diffusing particle, 
the diffusion will be confined to an interval of a finite length. Then the survival probability decays exponentially 
and, consequently, the mean exit time becomes finite[27|. Consider now our SFD model with N = 2. Then, for the 
left particle, the right one plays the role of a moving reflection boundary. That is, the left particle diffuses within 
a finite interval of a fluctuating length. Contrary to the case of the static boundary, the survival probability of the 
left particle decays as 1/i, and its mean exit time is infinite. It is interesting that the case iV = 3 brings us to still 
another behavior of the leftmost particle. For this particle, the central particle again represents a moving boundary. 
This boundary itself diffuses and, moreover, it feels another moving boundary to the right of it. As a result of such 
right-hand confinement the escape process of the leftmost particle is accelerated, its mean exit time is finite, but the 
variance of the exit time still diverges. For a general A^, the exit time of any n-th tagged particle with n > 4 possesses 
both a finite mean value and a finite variance. 

In leaving the subject, notice that the limit k — > oo alters only the prefactor A but not the long-time decay exponent 
of Sn{t). This limit corresponds to perfect absorption (particle is absorbed immediately when it hits the absorbing 
boundary). Consequently, lim 0„(t) gives the distribution of the n-th particle's first-passage time and lim (i„) 

equals the mean first-passage time. 

VI. CONCLUSION 

This paper provides a simple and general formalism for ID SFD problems with absorbing boundaries. Our main 
message concerns the relation of the physically motivated SFD problem and a broadly developed probabilistic field of 
order statistics and it can be summarized as follows. If the initial condition in the physical SFD problems is assumed 
in the form of a joint PDF for the order statistics based on the system of identical and independent parent random 
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variables, then this form is preserved during the subsequent time evohition. Implicitly, the observation has been 
already contained in the previous theoretical SFD studies in problems without boundaries. It can be easily read from 
the permanent-like structure of the Green function. 

Our work elaborates the connection between the SFD problem and the order-statistics point of view into new 
directions. First of all, we have proven that the relation continues to be valid in the problems with boundaries. In 
these problems, the arbitrary-time joint PDF for the positions of the particles has again the form of the joint PDF for 
the order statistics, provided one simply admits that the PDF for a single parent random variable has an atom at the 
origin. Secondly, specifically in the problems with boundaries, the joint PDF for the order statistics at an arbitrary 
time includes, in a mixed form, both the information on the particles remaining still in the pore and on their exit 
times. We have shown that these two pieces of information can be separated thereby yielding new exact results on 
the tagged particle diffusion and on the escape process, respectively. 

The robustness of the relation is easily understood as it lies in both cases (with and without the boundaries) in 
the hard-core inter-particle interaction. It is therefore not surprising that our results are independent on a detailed 
form of the single-particle dynamics. They hold true, for example, if the single particle is acted upon by a general 
space-dependent or even a space-and-time-dependent external force. The overall structure of our solution remains 
unchanged if we only change the details of the escape process and switch, for example, from the radiating boundary 
to the totally absorbing one. Physically motivated calculations in these directions are in progress and will be reported 
elsewhere. 

Eventually, let us comment on our choice of the initial condition as given by Eq. (|17p . Obviously, it is not the 
most general option. For example, what happens if the particles are initially placed at specific set of coordinates 
2/1 < 2/2 < ■ • ■ < J/TV? Taking a general initial condition, the partial PDF as defined by Eq. (PT|) does not reduce to the 
simple product form ((22)) . Instead, it assumes a rather involved permanent-like form similar to the RHS of Eq. (|2ip . 
With the necessary changes, one can carry out all the steps which give the partial PDFs /(^"'^'(xTv-fcli), the tagged 
particle PDF ^/'„(a;;t), and the survival probabilities Sn{t). But again, all the final expressions contain permanents. 
In spite of these more complicated finite-time expressions, the long-time dynamical exponents coincide with those 
derived in the preceding sections. Again, the parameters which specify the initial condition (together with D and k) 
appear in the prefactors of the power laws (03]), PSJ) . 
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Appendix A: Integration over the phase space TZn 

Let G(xi,X2, ... ,xn) be a symmetric function of its arguments, i.e., let 

G{xi,X2,.-,xn) = G(a;^^^(i),x^^^2),---,a;^„(Ar)) 

holds for any of Nl permutations 7r,„ of the subscripts {1, ..., N}. Then the integration of G{xn) with respect to 
the coordinates x^ over the phase space TZn (defined in Sec. IIII Cp extends to the integration with respect to the 
components of xjv over the half-lines as follows, 

dxNG{xi,X2,--- ,xn) — / da;i/ dx2 ... / dx^ G{xi,X2, ■■■ ,xm) (A1) 

Jo Jo Jo 

"R-N 

1 /• + 00 (■+00 /• + 00 

^ ~m j '^^^ I ^^2--- / dxN G{xi,x2, ... , xn) . 

The first equality introduces just another notation for the integration over TZn, the second equality follows from the 
symmetry of G{xn)- For a further discussion see, e.g.. Appendix C in Ref. l2ll 

In particular, if G{xn) stands for a product of identical functions, i.e., G{xn) = g{xi)g{x2)... g{xN), then we have 



AxNg{xi)g{x2)...g{xN) = -^ 



TZt. 



N 

Axg{x) 
Uo 



+ 00 

(A2) 
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If we are interested in properties of a tagged particle, say of the n-th one, we should carry out another type of the 
integration over the phase space TZn- In this case we have to integrate over the variables Xi, i ~ I,---, n—l,n+l,..., N. 
Let us denote the domain of such integration as 7?.Ar\x„. Further, let us define the symbol 

[dxN\xn] = da;i...da;„_ida;„+i...da;Ar . (A3) 

The integration of G{xn) = g{xi)g{x2)--- g{xN) over TZpf\xn runs as follows 

[<lxN\xn] 9{xi)g{x2)-.-g{xN)= g{xn) dxi g{xi) dx2gix2)... / dxn^ig{x„^i) X 



rXn rXn + i rXM-i 

X dXn+ig{Xn+l) / dXn+2g{Xn+2) ■■■ / dXNg{xN) 
Jo JO Jo 

r+oo 

dxg{x) 



g{xn) ' '•"" ^ "^^ ' '•+°° ' '''' 



NlINr 



dxg{x) 





(A4) 



where Nl = N — n {Nn = n — 1) stands for the number of particles located left (right) from the n-th particle (cf. 
Fig. [Ij. For a more comprehensive discussion of this integration technique, see Sec. Ill B. in Ref. USk 

Appendix B: Convolution of GF 

First, notice that the sum on the RHS of Eq. (|2T|) can be rearranged as 

p'^^\xN]t\yN]t') ^ ^p{x^^(^iy,t\yi;t')p{x^^(^2);t\y2;t')...p{x^^(^Ny,t\yN;t') . (Bl) 



Further, notice that p^^'{xN\ t \ y^; 0) is symmetric under the permutations of ijN- Inserting the initial condition ()17p 
into Eq. ((2T|l gives 



N 

f^^HxN;t)^mJdmp^^\xN;t\m;0)Ylfiyn). 



17 "=1 



We now insert p^^'{xN',t\yN',0) as given by Eq. (jBip . interchange the summation and integration, and use the 
symmetry of p(^)(xAr;i I 27jv;0) (cf. Eq. ((JT|) ). We obtain 



/(^^(fAr;t)= V dyN p{x^^^Ny,t\yN;0)fiyN) 

X / dyAr„i p{x^^(^N-i);t\yN-i;0)f{yN-i)--- 
Jo 

X dyi p{x^^^(^iy,t\yi;0)f{yi) 
Jo 

N\ 

= ^P{^■^rr^(N)^,t)p{x^^(N-l)■,t)...p{XJr,„(l)■,t) , 
771—1 

where the last equality follows from the definition (|4]). The last sum comprises A^! identical summands and hence we 
arrive at Eq. (P^ . 

Appendix C: Derivation of Sn{t) 

Let us define the auxiliary function B{x,t) = S{t) — F{x,t) with the properties: B{0,t) — S{t), B{+oo,t) = 0, 
p{x]t) = —■^B{x,t). Then we can rewrite Eq. ([55)1 as 

Sn{t) = -nf^\ £dx ^B-\l Bf- . (CI) 
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Using the binomial theorem we expand (1 — B)^~^ and we insert the unity in the form (n + k)/{n + k). Then we get 

^"W--L EViT . /d-^r-, (C2) 



where the definite integral equals {S{t)) . On the whole, we have proven Eq. (15^]) from the main text. Moreover, 
the sum in Eq. ([M)) can be represented by the definite integral (Sq 

Suit) =n[ ] I dz z"-i(l - z)^-" . (C3) 



n 



The integration by parts directly yields Eq. (|40j 
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